The DN, 7rS c interaction in finite volume and the A c (2595) resonance 
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In this work the interaction of the coupled channels DN and 7rE c in an SU(4) extrapolation of 
the chiral unitary theory, where the A c (2595) resonance appears as dynamically generated from 
that interaction, is extended to produce results in finite volume. Energy levels in the finite box 
are evaluated and, assuming that they would correspond to lattice results, the inverse problem of 
determining the phase shifts in the infinite volume from the lattice results is solved. We observe 
that it is possible to obtain accurate 7r£ c phase shifts and the position of the A c (2595) resonance, 
but it requires the explicit consideration of the two coupled channels. We also observe that some of 
the energy levels in the box are attached to the closed DN channel, such that their use to induce 
the 7r£ c phase shifts via Liischer's formula leads to incorrect results. 

PACS numbers: 11.80.Gw, 12.38.Gc, 12.39.Fe 



I. INTRODUCTION 

One of the topics where efforts are recently devoted 
within Lattice QCD is the determination of hadron spec- 
tra, both in the meson and baryon sector [H-[2l|. Af- 
ter earlier claims of a successful determination of the 
hadron spectra using rough approximations and large 
pion masses, work continues along this line with more ac- 
curate approaches and problems are arising that were not 
envisaged at first glance. The "avoided level crossing" is 
usually taken as a signal of a resonance, but this criteria 
has been shown insufficient for resonances with a large 
width [H-il]. The use of Liischer's approach [Mill is 
gradually catching up. It is suited for the case when one 
has resonances with one decay channel in order to pro- 
duce phase shifts for this decay channel from the discrete 
energy levels in the box. Yet, most of the hadronic res- 
onances have two or more decay channels and the need 
to go beyond Liischer's approach becomes obvious. This 
method has been recently improved in Ref. [24[ by keep- 
ing the full relativistic two body propagator (Liischer's 
approach has an exact imaginary part but makes approx- 
imations on the real part) and extending the method to 
two, or more coupled channels, which had also been ad- 
dressed before |27H29| . and continues to catch the atten- 
tion of the practitioners [3(1 Hl|. The new method is 
also conceptually and technically simpler and serves as a 
guideline for future lattice calculations. Continuation of 
this new practical method have been done in Ref. (32j for 
the application of the Jiilich approach to meson baryon 
interaction and in Ref. [33| for the interaction of the DK 
and r]D s system where the D s *o(2317) resonance is dy- 
namically generated from the interaction of these par- 
ticles [34l - l37| . The case of the k resonance in the Kir 
channel is also addressed along the lines of Ref. [24| in 



Ref. [H. 

The investigation of different problems following 
Refs. [3, HH, HH is showing that every case stud- 
ied reveals particular features and there is no common 
behaviour in the levels of the finite box nor on the way 
that the phase shifts or bound states are obtained from 
the spectra of the finite box. The study of these problems 
along those lines is most useful, since it sheds light on how 
to deal with results of QCD lattice calculations, which 
precision is needed in the lattice results to accomplish 
a desired accuracy in the phase shifts or resonance pole 
positions, and which strategy is most useful to gather 
lattice results from where the results in infinite volume 
can be obtained with maximum accuracy. 

The case we report here is one more in this line, show- 
ing as a novel feature how in this case some of the low 
lying levels in the finite box, for energies where the 7r£ c 
channel is open and the DN one is closed, are tied ba- 
sically to the DN channel, such that the blind use of 
Liischer's approach would lead to unrealistic phase shifts 
for the 7rE c channel and the resonance position. We show 
the problems that arise in this analysis and provide the 
two channel approach to solve them. These results should 
be most useful when QCD lattice results are produced to 
describe the A c (2595) resonance. 

This paper is organized as follows. In Sec. II we show 
the formalism of the DN and 7rE c interaction in infinite 
and finite volume. In Sec. Ill, the inverse problem of 
getting the phase shifts from two channels analysis is 
shown, while in Sec. IV the results obtained by using one 
channel analysis are shown. Finally, a short summary is 
given in Sec. V. 



II. FORMALISM 



A. The DN, 7rE c interaction in infinite volume 
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The DN system, in collaboration with coupled chan- 
nels, leads to the formation of a meson-baryon composite 
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state, the A c (2595) [39t - 442l | . In this section, we briefly re- 
visit the form of the 7rE c and DN interactions from the 
chiral unitary approach. This will allow us to review the 
general procedure of calculating meson-baryon scattering 
amplitudes. In the chiral unitary approach the scattering 
matrix in coupled channels is given by 



T = [l-VG\- l V 



(1) 



where V is the matrix for the transition potentials be- 
tween the channels and G, a diagonal matrix, is the loop 
function for intermediate 7r£ c and DN states, which is 
defined as 



G 



d 4 q 



2M 



(2tt) 4 q 2 



+ ie(P- q) 2 - M 2 + ie 



(2) 



where m and M are the masses of the tt or D meson and 
the baryon E c or N, respectively. In the above equation, 
P is the total incident momentum of the external meson- 
baryon system. 

In the present problem we have two main channels, 
7rE c and DN. There are also other channels considered 
in Refs. (40| - |42T| . such as rj A c , KE C , KE' C , D S A, and 77' A c , 



which play a negligible role to generate dynamically the 
A c (2595) state, and are not considered here. 

We study only the s-wave interaction, hence, the tran- 
sition interaction (potential) for channel i to j reads, 



Vu = 



.9^(2^- Mi -M S 



2Mi 



1/2 



2M, 



1/2 



(3) 



where / = 93 MeV is the pion decay constant, Ei and 
Ej are the energy of incoming/outgoing baryon E c or N. 
The transition coefficients Cij are symmetric with respect 
to the indices, and also isospin-dependent. By naming 
the channels, 1 for 7rE c and 2 for DN, the coefficients 
for the case of isospin I = are 



Cix — 4, C12 — C21 — — — , 



Cv> — 3. 



(4) 



The loop function G can be regularized both with a 
cutoff prescription or with dimensional regularization in 
terms of a subtraction constant. Here we make use of the 
dimensional regularization scheme. The expression for G 
is then 



Gf(s,m i ,M i ) 



+ 



2Mj 
(4tt) 2 

Q*(v^) 



772 - 



2.s 



, Mr 



(5) 



[log (a - {Ml - m 2 ) + 2^Q,(V^)) + log (a + (Mf - mf) + 2^Qi(^)) 



log (-a + (Mf - ml) + 2^- s Q l (V7s)) - log (-s - (Mf - mf) + 2 v / iQ l (v^))] 



where s = E 2 , with E the energy of the system in the 
center of mass frame, Qi the on shell momentum of the 
particles in the channel, /i a regularization scale and cij (fi) 
a subtraction constant. The form of Eq. Q is adapted 
from the light hadron sector to the charm sector using 
SU(4) symmetry but reducing the strength of the dia- 
grams with a heavy vector exchange by the weight of 
its propagator. One should, in principle, not expect a 
good SU(4) symmetry, and actually it is broken in this 
case through the use of physical masses in the propa- 
gator. But in basic vertices the symmetry works quite 
well (see an extensive review in Ref. [43[ and references 
therein, and in section IV of the paper E3). We thus 
follow this approach, as done in Refs. [39h42| | . One must 
bear in mind that uncertainties in the actual value of 
Vij in Eq. (J3]) are taken into account, in the spirit of the 
renormalization group, by means of the freedom in the 
subtraction constants of G in Eq. ([5]), which are adjusted 



such as to get the energy of the A c (2595) in the right 

place. Note that the only parameter-dependent part of 

2 

G is a(/i) + ln^T-. Any change in /1 is reabsorbed by a 

change in a(/j) through a(/x ) — a(/i) = ln'V- 

In the infinite volume case, the use of Eq. ([IJ with 
the two channels that we consider leads to a dynami- 
cally generated state at the energy of 2596 MeV, which 
we associate to the A c (2595) resonance, using the dimen- 
sionally regularized G D function with fi — 1000 MeV 
and the subtraction constant a = —2.02, which are val- 
ues of natural size. In Fig. [T]the modulus squared of the 
DN — » DN scattering amplitude as a function of the in- 
variant mass of the DN system for Ion — is shown. 
There is a clear rather narrow peak around 2596 MeV 
which indicates the state of A c (2595) as a bound state of 
DN. 

The poles of the amplitudes are found in the second 
Riemann sheet, where in the loop function G D in Eq. (J5]) 
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FIG. 1: Modulus squared of the DN — > DN scattering am- 
plitude for Idn — 0. 



constant gfj N that we obtained above and the loop func- 
tion G% N , we find that g 2 DN ^eL = _o.91 - i0.08 at 
E = y/sRi and, thus, about 91% of the sum rule comes 
from the DN state, indicating that we have largely a 
bound DN channel in our approach. 

Another way to check how important is the DN chan- 
nel in generating dynamically of the A c (2595) state is to 
change the parameters of the potential slightly, making 
the A c (2595) disappear. This can be achieved, for exam- 
ple, by merely reducing the strength of the potential V22 
that describes the scattering in the DN — > DN chan- 
nel. Namely, we replace V22 T1V22 and vary 77 between 
1 and 0, and as can be seen in Fig. [2l the phase shifts 
for the 7rE c — > 7rE c scattering amplitude obtained from 
the Bethe-Salpeter equation, Eq. ([T]), drastically change 
with the strength of V22 ■ The normalization that we use 
is such that in one channel Hm 2 



Qi(yfs) is changed to —Qi(yfs) for the channel where 
Re-^/s is above threshold. The couplings of A c (2595) 
to the 7rE c and DN channels can be obtained from the 
residues at the pole, by matching the amplitudes to the 
expression 



T; 



9i9j 



(6) 



for y/s close to the pole yfsR. We find y/sR = (2596.3 — 
il.6) MeV. The couplings, g± and gj are complex in 
general. In this way, we get g 2 s = 0.65 — i0.10 and 
0^ = 43.3 + 43.6. 

In Ref. [45.], an analytical study of bound states in 
problems with coupled channels has been performed. A 
modern formulation of the compositeness condition of 
Weinberg (4(| for coupled channels has been derived from 
a sum rule that comes from the normalization to unity 
of the wave function of the bound state 1 , 



T(E) 



~8nE 



2ms c p cot 6(p) — ip 



(8) 



from where we determine the phase shifts in the infinite 
volume problem. 
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(7) 



E—E a 



where g 2 are the residues of the Tu scattering matrices 
(coupling squared) at the pole of the bound state (a) 
and Gu are the propagators of the two particles of the 
corresponding channels (the loop function G D that we 
use here). 

One can see from the derivation in Ref. [45[ that each 
term in Eq. (J7J accounts (with reversed sign) for the prob- 
ability of the bound state to be made by the pair of par- 
ticles of the channel considered. By taking the coupling 



FIG. 2: Phase shift for 7rE c scattering derived from the cou- 
pled channels unitary approach with different values of r\. 
Solid line: rj — 1. Dashed line: r\ = 0.8. Dotted line: only 
7r£ c channel considered. 

In Fig. [2 the solid curve stands for the phase shifts for 
7rE c scattering derived from the two coupled channels 
unitary approach with r\ = 1, while the dashed line is 
obtained with ij = 0.8. We also show with the dotted 
line the phase shifts that were obtained considering only 
the 7rE r channel. 



1 Technically, the state obtained is a resonance because the 7rS c 2 We mention that in the present calculation, we replace T(E) of 
channel is open, although with very small pase space. However, meson-meson system in Ref. by 2m£ c T(E) since we treat a 

the sum rule also holds at the resonant pole as shown in Ref. [47ll . meson-baryon system. 
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B. The DN, 7r£ c interaction in finite volume 

One can also use regularization with a cut off in three 
momentum once the q° integration is analytically per- 
formed 14a with the result 



performing an average for different q max values between 
2000 MeV and 4000 MeV one obtains a perfect conver- 
gence, as one can see in Fig. [3] Note that the imaginary 
part of G D and that of the integral in Eq. (TTTT) are iden- 
tical and they cancel in the construction of G, which is a 
real function. 



I«1<-3n 



d 3 q 



2m 2 



uii(q) + u> 2 (q) 
E 2 -(a Jl (q)+UJ 2 (q)) 2 + te 



(9) 



with mi, m 2 corresponding to m and M of Eq. ([5]) for 
each channel j. In Ref. (49| the equivalence of the two 
methods was established. 

When one wants to obtain the energy levels in the box, 
one replaces the G function by G, where instead of inte- 
grating over the energy states of the infinite volume, with 
q = | q | being a continuous variable, as in Eq. (J9j) , one 
sums over the discrete momenta allowed in a finite box of 
side L with periodic boundary conditions. We then have 
G = diag(Gi,G 2 ), where 



Gi = 



1 

L 3 

T 



l9l<9n 



E 



2m 2 



W\{q) + uj 2 (q) 



2u 1 {q)u 2 {q) E 2 — (ui(q) + w 2 (g)) 2 



n, n e Z , 



(10) 



with the same notation as in Eq. (J9)). 

By using the dimensional regularization of the loop 
function G D of Eq. ([5]), we can write [33| 



G(E) = G D (E) 



+ lim 

Qmax - >00 



[ L 3 

= G D (E) + B, 
where I(q) is the integrand of Eq. (JSJ 



q<q a 



(27T) 



2m 2 



2wi(g) w 2 (g) E 2 - (cj^q) + u 2 (q)) 2 + le 



(11) 



.(12) 



The three dimensional sum in Eq. (|llj) can be re- 
duced to one dimension considering the multiplicities of 
the cases having the same n 2 [32|, • The integral in 
Eq. (|11|) has an analytical form as shown in the appendix 
of Ref. [ll| (see erratum). 

In the box, the same Bethe-Salpeter equation is used 
substituting G D by G of Eq. (fTTj) . When calculating the 
limit of q max going to infinity in Eq. (jlll) we obtain os- 
cillations which gradually vanish as q max goes to infinity. 
Yet it is unnecessary to go to large values of q m ax, and 
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FIG. 3: Real part of the second part of Eq. ()lip for the 7rS c 
channel. The dashed line stands for the average that we take 
between 2000 MeV and 4000 MeV for q max . The results cor- 
respond to a value of L — 2.5 m^ 1 and E= 2600 MeV. 



The eigenenergies of the box correspond to energies 
that produce poles in the T matrix. Thus we search for 
these energies by looking for zeros of the determinant of 
1 - VG 



det(l - VG) 



1 — V\\G\ — V 22 G 2 

r2 



+{Vi 1 V 22 -V{ 2 )G l G 2 = 0. (13) 



In Fig. @] we show the first five energy levels obtained 
for the box for different values of L. We observe a smooth 
behavior of the levels as a function of L. The second level 
is special because it mainly comes from the DN channel. 
If we only include the DN channel, we get very similar 
results as shown with the dashed curve. From Fig. [4] we 
can also see that the second level is rather independent 
on the values of the cubic box size L. The value for the 
eigenenergy of the second level is around 2600 MeV which 
is very close to the mass of A c (2595). 



III. THE INVERSE PROBLEM OF GETTING 
PHASE SHIFTS FROM LATTICE DATA 

In this section we face the problem of getting bound 
states and phase shifts in the infinite volume from the 
energy levels obtained in the box using the two channel 
approach of Ref. , which we would consider as "syn- 
thetic" lattice data. To accomplish this we need more 
information than just the lowest level, but we shall see 
that the first two levels shown in Fig. 2] already provide 
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FIG. 4: Energy levels as functions of the cubic box size 
Z^m^ 1 ), derived from the chiral unitary approach. 
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FIG. 5: Energy levels as functions of the cubic box size L, 
reconstructed from fits to the "data" of fig. [4] by using the 
potential of Eq. (|14|l . The band corresponds to the statistical 
errors of the fitted parameters. 



the necessary information to reproduce the problem in 
the infinite volume. 

In Ref. (24[ several methods were suggested to solve 
the inverse problem, one of the methods was the follow- 
ing: In two channels one has three degrees of freedom, 
Vii, V12, and V22, or in terms of phase shifts, Si, S2 and 
the inelasticity 77 (or equivalently the mixing angle) . One 
strategy to obtain these magnitudes is to use three lev- 
els that contain a certain energy and using Eq. (13) 
determine the three degrees of freedom for a given en- 
ergy. This strategy is used in [20] and is also suggested 
in [13, HH to obtain directly Si, 82 and r\. The technical 
problem that this method poses is that for a given en- 
ergy one might need a too large or a too small value of L 
that could make the computation too lengthy or inaccu- 
rate, respectively, although using moving frames, like in 
[30, [31j one obtains more levels that can reduce the span 
of values of L. Other different methods were suggested in 
Ref. (24| and we borrow here the one based on a fit to the 
data in terms of a potential parameterized as a function 
of the energy suggested by the coupled channel unitary 
approach of the work of Ref. (3tJ or Refs. [33 - I36I ]. As we 
can see in Eq. (J3J, the potentials have a large constant 
part, some terms proportional to s and some terms pro- 
portional to y/s. It is very easy to see that if one chooses 
a region of energies around a certain value of s, So, the 
potential can be expanded as a function of s — sq to a 
good approximation. Choosing sq — (m, + M^ c ) 2 then 
the ansatz of the following equation 

Vij = o« + bij (s - (m, + M Sc ) 2 ) (14) 

is a very accurate assumption. 

We assume that the lattice studies provide us with 
ten eigenenergies corresponding to the first two levels 
of Fig. 0]for different values of L between 1.7 m^ 1 and 
3.3 m" 1 . We also assume that the levels are provided 
with an error of ±10 MeV. We make a best fit to the 



data assuming a potential as in Eq. (fill) . We look for 
the minimum \ 2 by the MINUIT fit program and obtain 
a set of parameters for , bij . The minimum % 2 that 
we get from the best fit is 2.5 x 10~ 2 . Then we generate 
random sets of the parameters within the range of error 
of each parameter determined by the best fit, such that 
X 2 is only increased below Xmin + 1- With these values 
we generate the spectrum of Fig. |4] by searching for the 
zeros of the determinant of 1 — VG. This provides a band 
of values for the spectrum shown in Fig. [5] As found in 
Refs. [H, [H, [H, HI] one has the freedom to choose the 
regularization constant different to the one used origi- 
nally to generate the spectrum and the good fit to the 
spectrum is obtained by a corresponding change in the 
parameters of the potential, a feature tied to the renor- 
malization group. This is most welcome because, when 
the lattice data are provided to us, we do not know which 
implicit regularization subtraction constant the lattice 
data is supporting (the lattice spacing is not a problem in 
this sense, as discussed in Ref. [32| )■ The inverse method 
has only a real value if the results that one obtains are 
independent of this subtraction constant. 

The aim of the inverse method is to get the phase shifts 
and bound state in the infinite volume from the spectrum 
obtained in the box. For this purpose we take now the 
potential obtained by the best fit to the synthetic data 
with a chosen subtraction constant, and use it in Eq. ([1]) 
to produce the scattering amplitude in the infinite volume 
case using G D with the same subtraction constant. 

In Fig. [6]we show the 7rS c — > 7rS c phase shifts, for I = 
0, Sq, obtained with the uncertainty provided by the set 
of parameters that fulfill the x 2 < Xmin + 1 condition. As 
we can see, the agreement with the exact results is quite 
good, and we see how the errors in the determination of 
the lattice levels have propagated in the determination of 
the phase shifts. The results obtained show the presence 
of the A c (2595) resonance around this energy with an 



6 



uncertainty of ±5 MeV. The width that we calculate from 
the position of the pole in the complex plane is about 3 
MeV. One can also obtain this from the phase shifts by 
using 



T = 



9 



1 



-8ttE 



E - E R + iT/2 2m Sc pcot8(p) - ip 
which leads to the equations 

E - E R -ms c p cot S (p) 



9 A 



4nE 



r 



AttE 



(15) 

(16) 
(17) 



from where we get T = 3.3 ± 1.9 MeV. 



300 




2600 2620 2640 2660 2680 2700 
E(MeV) 

FIG. 6: Phase shift for 7rE c scattering obtained from the two 
coupled channels unitary approach (solid line). The band is 
derived from the fits to the "data" of Fig. [5] by using the 
potential of Eq. (|14[) with two channels. 
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FIG. 7: As in Fig. O but for the new fit with the error ±5 
MeV for the centroids of the 'data'. 



bij{y/s- (m„ 



M Sc ) 2 ), 



(18) 



hence, introducing three more parameters. We should 
note that with such large number of parameters one will 
not only be producing a smooth fit to the data, but the fit 
will also search situations to reproduce the fluctuations. 
This means we should consider the extra errors that come 
from this source certainly as an upper bound. The new 
results, considering also the dispersion of the centroids, 
are shown in Fig. [5] (note that we added also a new data 
point). As we can see, and similarly to what was fond in 
Ref. [3g], the error band has increased and now extends 
over the whole range of the assumed errors of the lattice 
data. 



At this point we want to improve the error analysis 
by introducing two new ingredients. The first one is to 
consider that the centroids of the data are not exactly on 
the exact curve, as it would correspond to actual lattice 
data. For this purpose we follow exactly the procedure 
done in Ref. [2J] and let the centroids move randomly 
within ±5 MeV from the exact point in the curves of 
Fig. fJJ We then make a large number of runs and deter- 
mine the new band of results. This is shown in Fig. [7] 
As we can see, the error band has increased a bit, by 
about 20%. We repeat the procedure allowing now the 
centroids moving randomly with ±10 MeV and we find 
that the band increases by an extra 20% with respect to 

Fig. m 

In a second step we also want to take into account the 
effects of using more freedom in the parameterization of 
the potential. This was discussed in Ref. [38[ and was 
shown to be an extra source of uncertainty. We want to 
investigate what happens here. For this purpose we use a 
different parameterization of the potential guided by the 
terms that appear in Eq. Thus, we take 
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As in Fig. [H but for the new fit by using the potential 
l8l) with two channels. 



FIG. i 
of Eq. 



With the new potential we reevaluate the phase shifts 
with two channels as we have done for the Fig. [51 and the 
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new results are shown in Fig. [9] The agreement with the 
exact results is good but the error band is now increased. 
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FIG. 9: As in Fig. [6j but the band is derived from the fits to 
the "data" of Fig. [8] by using the potential of Eq. (JTHJl with 
two channels. 



are below the 7rS c threshold. Besides, the second level is 
very stable and only gives us an energy point, with errors 
in the energy bigger than the width of the resonance. It 
is not possible to reconstruct the phase shifts in these 
circumstances. Because of this, it is more appropriate to 
consider all the data of one level, since one is then using 
the information on the correlation of these data. Then 
we fit all these data (five points) with the potential of 
Eq. (fT4|) but with only the 7rS c channel (2 parameters). 
We get a good fit to the first level with Xmin = 3.4 x 
1CP 2 . Then we use this potential to determine the infinite 
volume scattering amplitude T 



T(E) 



1 



V- 1 - G{E) ' 



(21) 



The results for 7r£ c phase shifts, with the lower band, 
are shown in Fig. [TO] As we can see, the results obtained 
differ substantially from the exact results obtained with 
the two coupled channels, as shown in Fig.[10]by the solid 
line. We should note that, not even the scattering length 
is provided correctly. 



The analysis done here is based on the dominance of 
the 7rS c and DN two body channels. The inclusion of 
three body channels makes the work technically much 
more involved, but has already been addressed formally 
in Ref. [Hj]. In the present case the A c (2595) resonance 
can decay into 7rS c but also in uncorrelated A c 7T7t. For- 
tunately the branching ratio in this channel is small and 
with large uncertainties, 18 ± 10% [53j], that justifies its 
neglect in the present work. Should there be better data 
in the future on this branching ratio, and should one aim 
at a very accurate solution, some work along the lines of 
Ref. [53 would be advisable. 



IV. ONE CHANNEL ANALYSIS 

Since the 7rE c channel is the only one open, one might 
be tempted to apply Liischer's approach with just one 
channel. One would be assuming implicitly that the ef- 
fect of the DN channel would be absorbed in the 7rE c 
potential. In such a case, the energy spectrum of Fig. [5] 
would be given by the poles of 



T(E) = 



1 



V- 1 - G{E) 



(19) 



which gives us V^ 1 — G for a value of E eigenenergy of 
the box. For this energy, we can then write the scattering 
amplitude in infinite volume as 



T(E) 



1 



G(E) - G(E) 



(20) 



However, the direct application of this formula does 
not give information above the 7rE c threshold from the 
first level, since the eigenenergy values of the first level 
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FIG. 10: Phase shift 6® for 7rE c scattering. Solid line: exact 
results obtained from the two coupled channels unitary ap- 
proach. Lower band: derived from the fits to the first level 
"data" of Fig. [5] by using the potential of Eq. (|14|) with only 
one 7r£ c channel. Upper band: derived from the fits to the 
second level " data" of Fig. [5]by using the potential of Eq. (|14|) 
with only one 7rE c channel. 



Similarly we take now the second level in Fig. [5] and 
perform the same exercise. The best fit gives a larger 
Xmin = 0.4. The phase shifts obtained are also shown in 
Fig. [TU] with the upper band. They are also in very bad 
agreement with the real results. 

Besides, we can take the second level of Fig. [5] and 
analyze it with just the DN channel. Then we cannot 
get the 7rE c phase shifts, but we see that the scatter- 
ing amplitude obtained in the infinite volume by using 
Eq. ([21} has a pole at E = 2601 ± 6 MeV, correspond- 
ing to a DN bound state. This is telling us that the 
second level of Fig. [5] is mostly tied to the DN channel. 
One could guess that the stability of the energy level 
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as a function of L is showing the presence of a bound 
state, although this is not always the shown in 

Ref. [3J|. Yet, in the present case we have gone one step 
beyond, because our inverse analysis provides a param- 
eter set for the potential from where we can determine 
poles, couplings and then test Eq. ([7]) for the composite- 
ness condition. In the present case this renders a value of 
around 0.9 for —g DN dG d '^ N , from where we conclude that 
the state corresponds essentially to a DN bound state, 
weakly decaying into 7rS c . 

Once again we repeat the exercise done at the end of 
the former section and take into account at the same 
time the two effects considered here. Thus, we consider 
the dispersion of the centriods of the data and use the 
new parametrization of Eq. (|18|) . In this case we have 
now three parameters in the potential. The new results 
can be seen in Fig. [TTJ We can see that the quality of 
the fit to the two levels is worse than the one obtained in 
Fig. [5] with two channels. 
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FIG. 12: As in Fig. [10] but, the band is derived from the fits 
to the first two levels "data" of Fig. [11] by using the potential 
of Eq. (|18[) with only one 7rE c channel. 
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FIG. 11: As in Fig. [5] but for the new fit with the error ±5 
MeV for the centroids of the 'data' and by using the potential 
of Eq. (|18[) with only one 7rE c channel. 



Next, we evaluate the phase shifts with this new po- 
tential, fitted to the two levels, and the results can be 
seen in Fig. 1121 Once again, we see that assuming just 
the 7rE c channel in the analysis leads to unrealistic 7rS c 
phase shifts. 

In summary, we see the problems that arise when we 
try to use Luscher's approach for the interpretation of 
the lattice spectrum in a case where the relevance of a 
closed channel is huge, like in the present case. We also 
see that doing the fit analysis, that allows us to circum- 
vent Luscher's approach, but using only one channel also 
fails to provide realistic phase shifts. The analysis with 
two channels shows here the tremendous power that the 
coupled channel approach has in this case. 



V. A TEST IN TERMS OF A CDD POTENTIAL 

On the other hand, there is the possibility that the 
nature of A c (2595) resonance could be a genuine state, 
not dynamically generated by the DN interaction. For 
this purpose we have made a test introducing a differ- 
ent potential for the DN interaction where a CDD pole 
(Castillejo, Dalitz, Dyson) 54] is introduced by hand. 
The potential for DN interaction now is, 



V = V M + r yCD ° , (22) 

V s - \/SCDD 

where Vm is assumed to be energy independent and 
9c ddi scdd are the parameters of the CDD pole. 

Like it has been done in Ref. (33[, from the above po- 
tential, we now find 

1 _ 9cdd g2 = 1 (93) 

U/i-V^F) 2 f i-z' 

with Z the field renormalization constant for the gen- 
uine state, which accounts for the probability to have a 
genuine state. 

We take the potential Eq. ([22]) with Vm of the or- 
der of 10 times smaller that the potential used for V22, 
\/scdd corresponding to a 20 MeV below the mass of 
the A c (2595) and then g% DD = 1.86 such as to get the 
bound state at the mass of A c (2595). We find that at 
the pole of this state, from Eq. (|2"5]l . we get Z = 0.96, 
which shows that the introduction of a CDD pole as in 
in Eq. (f2"2"f is good enough to generate a genuine state. 

The energy levels in the box with the potential of 
Eq. (|2"2")1 for the DN channel are shown in Fig. [T3] As we 
can see, the levels are different from those obtained with 
the couple channel potential, which are shown in Fig. 2] 
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It is clear that the determination of the levels with lattice 
calculations can differentiate between the two different 
types for the potentials. 



3100 r 
3000 ^ 
2900 ; 

> 

2800 - 

LU 

2700 : 
2600 '-- 



2500 r 1 

1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 

FIG. 13: Energy levels as functions of the cubic box size 
L(m^ 1 ), derived from the potential of Eq. (1221) for the DN 
channel. 

Next, we perform a fit, with a CDD potential Eq. (j2"2")l 
for V22, to the first two levels of Fig. |4] that are obtained 
with the couple channel potential. We have taken also 
ten points over the curves and assumed ±10 MeV er- 
rors as we have done above. The best fitting results are: 
9c d d = 91.4 ± 26.9, y/Ic~DD = 8190 ± 1600 MeV, and 
Vm = 1.23 ± 0.45, with these values and their uncertain- 
ties, we can get now, at E = 2596 MeV, Z = 0.15 ± 0.07 
from Eq. (|2U1) . which is in the order of 15% with large er- 
ror. This means around 85% fraction of A c (2595) being 
dynamically generated. We should note that although 
we get a good fit with a potential that formally contains 
a CDD pole, the large value of the mass of the CDD 
pole renders the potential smooth, as in the case of the 
coupled channels, and the test tells us that the state cor- 
responds to a dynamically generated one. We should 
also mention that we now only take two energy levels of 
Fig. @] for fitting, this is why the best results have large 
errors. If we took more energy levels, we could determine 
these values with more precision. But our present result, 
Z = 0.15±0.07, suffices to show that from this limited 
information one can get valuable conclusions on the na- 
ture of the A c (2595) resonance, which is mostly a DN 
bound state. 



VI. SUMMARY 

In this work, we study the interaction of the coupled 
channels DN and 7rE c in an SU(4) extrapolation of the 
chiral unitary theory. The resulting interaction is used 
to reproduce the position of the A c (2595) resonance in 



the isospin zero DN channel. Then we conclude that the 
A c (2595) is mostly a DN bound state. 

We then study the interaction of the coupled channels 
DN and 7rE c in the finite volume. Energy levels in the 
finite box are evaluated. We assume that the results ob- 
tained would correspond to results given by lattice cal- 
culations. From there we address the inverse problem. 
We propose a rather general and realistic potential and, 
using two coupled channels, a fit to the synthetic data is 
made assuming some reasonable errors in the data. Then 
this potential is used in the infinite volume case, gener- 
ating the 7rE c phase shifts within an error band around 
the original results. This part provides information for 
lattice QCD calculations about the accuracy in the ener- 
gies of the spectrum needed to get a desired accuracy in 
the phase shifts. 

A second part of the investigation was about the use 
of a one channel Liischer's approach, with just the open 
7rE c channel, to induce 7rE c phase shifts from the finite 
volume spectrum. We found in this case that, due to 
the large weight of the closed DN channel in this prob- 
lem, the results obtained using Liischer's approach with 
just the 7rE c channel was of no use. Even more, mak- 
ing a fit analysis to the lattice data with just the 7rE c 
channel produced erroneous 7rE c phase shifts. Certainly 
one does not know a priori from the lattice QCD results 
whether two channels would be necessary in the analysis. 
However, we also showed that the results obtained from 
the analysis of the first two levels with just one chan- 
nel were different to each other. This could be taken as 
a clear indication that at least two channels are needed 
in a realistic analysis of the lattice QCD results in such 
a case. The results from the chiral unitary approach, 
and the two channel formalism shown here, which can 
be trivially generalized to more channels, provide a good 
perspective to undertake future lattice QCD calculation 
in this sector. We also showed that the analysis done 
here, not only provides us with the 7rE c phase shifts and 
the presence of a bound state, but through the test of the 
sum rule of Eq. (essentially Weinberg's compositeness 
test), it also tells us that this bound state corresponds to 
a molecular state of a DN system. 
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